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Abstract 

Using data taken with the CLEO III detector at the Cornell Electron Storage Ring, we have 
investigated the direct photon spectrum in the decays T(1S) ^ ygg, T(2S) ^ ygg, T(3S) ^ 'jgg. 
The latter two of these are hrst measurements. Our analysis procedures differ from previous 
ones in the following ways; a) background estimates (primarily from tt® decays) are based on 
isospin symmetry rather than a determination of the tt® spectrum, which permits measurement 
of the T(2S) and T(3S) direct photon spectra without explicit corrections for vr® backgrounds 
from, e.g., Xbj states, b) we estimate the branching fractions with a parametrized functional form 
(exponential) used for the background, c) we use the high-statistics sample of T(2S)^ 7r7rT(lS) 
to obtain a tagged-sample of T(1S)^ y + X events, for which there are no QED backgrounds. We 
determine values for the ratio of the inclusive direct photon decay rate to that of the dominant 
three-gluon decay T ^ ggg (R^ = B{ggy)/B{ggg)) to be; R^(1S)=(2.70±0.01±0.13±0.24)%, 
R..y(2S)=(3.18±0.04±0.22±0.41)%, and R^(3S)=(2.72±0.06±0.32±0.37)%, where the errors shown 
are statistical, systematic, and theoretical model-dependent, respectively. Given a value of Q^, one 
can estimate a value for the strong coupling constant as{Q‘^) from R^. 


*Submitted to the 33''*^ International Conference on High Energy Physics, July 26 - August 2, 2006, Moscow 
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Introduction 


Production oi a. BB meson pair, the Zweig-favored decay mode of T mesons, is not 
energetically possible for resonances below the T(4S), thus the decay of the T(1S) meson 
must proceed through Zweig-suppressed channels. Since the charge conjugation quantum 
number of the T resonances is C=~l, the three lowest-order hadronic decay modes of the 
T(1S) meson are those into three gluons (ggg), the vacuum polarization QED decay T ^ gg, 
and two gluons plus single photon (ggy). For the T(2S) and T(3S) resonances, direct 
radiative transitions, both electromagnetic and hadronic, compete with these annihilation 
modes. Since Tggg oc and Tgg^ oc a'^aem, the ratio of the decay rates from these two 
processes can be expressed in terms of the strong coupling constant [l|: 

^ + (2.2 ± 0.8)a.A]. (1) 

^ggg ^ggg ^ 

In this expression, the bottom cmark charge g?, = —1/3. Alternately, one can normalize to the 
well-measured dimuon channelj^ and cancel the electromagnetic vertex: {Tgg^/Tgg) cx af. 

In either case, one must dehne the value of appropriate for this process. Although 
the value ~ seems ‘natural’, the original prescription of Brodsky et a/.P] gave 

= (0.157Mx(is))^ for T(1S)^ ggy. 

Theory prescribes the differential spectrum (PN/dx^d cos 9z = p-y/Fbeam, and cos 6*2 

is defined as the polar angle relative to the e+e“ beam axis). The limited angular coverage 
of the high-resolution CLEO III photon detection, as well as the large backgrounds at low 
momentum due to decays of neutral hadrons (primarily 7 r°, rj, r]' and a;) to photons, plus 
the large number of radiated final state ‘fragmentation’ photons in this regime limit our 
sensitivity to the region dehned by jcos^^l <0.7 and x^ >0.4. We must therefore rely on 
models for comparison with the observed direct photon spectrum and extrapolation of the 
direct photon spectrum (excluding the fragmentation component) into lower-momentum and 
larger polar angle regions. 

Originally, the decay of the ground-state vector bb bottomonium into three vectors (both 
T —> ggg as well as T ^ ggi) was modeled in lowest-order QCD after similar QED decays of 
orthopositronium into three photons, leading to the expectation that the direct photon spec¬ 
trum should rise linearly with x^ to the kinematic limit {x^ Ih phase space considerations 
lead to a slight enhancement exactly at the kinematic limit j3|,|^. Koller and Walsh consid¬ 
ered the angular spectrum in detailj^, demonstrating that as the momentum of the most 
energetic primary parton (photon or gluon) in T —ygg or T —>■ ggg approaches the beam en¬ 
ergy, the event axis tends to align with the beam axis: x^ ^ 1 ^ dN/d{cos6z) 1 - 1 -cos^ 6^. 
Field 0 argued that = 1 is non-physical, since it corresponds to a recoil gg system with 
zero invariant mass, while the recoil system must have enough mass to produce on-shell final 
state hadrons. Using a phenomenological parton shower Monte Carlo technique which took 
into account the correlation of photon momentum with recoil hadronization phase space. 
Field predicted a significant softening of the lowest-order QCD predicted spectrum, with a 
photon momentum distribution peaking at x.y ~ 0.65 rather than x.y —> 1. (In the limit of 
completely independent fragmentation, the same argument, in principle, would apply to 3- 
gluon decays.) This result seemed in conflict with the extant CUSB0 data, which indicated 
a spectrum more similar to the lowest-order QCD prediction. A subsequent measurement 
by CLEO-lj3], however, favored Field’s softened spectrum over lowest-order QCD. Given 
the poor resolution of the CLEO-I electromagnetic calorimeter, that measurement was also 


3 





consistent with a snbseqnent modification to lowest-order QCD which calcnlated corrections 
at the endpoint by snmming leading logs of the form ln(l-a:^). Higher statistics measnre- 
ments by Crystal Ball0 as well as ARG US flil corroborated this softened photon spectrnmd 
A snbseqnent CLEO analysis (CLEO-II)|ll|, based on ~1 M T(1S) events, provided a high- 
statistics confirmation of a photon spectrnm peaking at ~ 0.65, and was able to trace 
the direct photon momentnm spectrnm down to ?«0.4; at that momentnm, the direct 
photon signal becomes less than 10 % relative to the backgronnd, whereas the systematic 
errors on the backgronnd estimate in that momentnm region exceed 10%.^ Contemporary 
with the CLEO-II analysis, Catani and Hantmann first pointed ont complications dne to 
the presence of fragmentation photons emitted from final-state light qnarks downstream of 
the initial heavy qnarkonia decay [l^ (essentially final state radiation). These can dominate 
the backgronnd-snbtracted spectra for x^ <0.4 and therefore (if not corrected for) lead to an 
over-estimate of the T( 1 S)^ gg''y branching fraction, and an nnderestimate of the extracted 
valne of 

Hoodbhoy and Ynsnfjl^ also performed a rigorons calcnlation of the expected T(1S)—> 
gg'j decay rate, by snmming all the diagrams contribnting to the direct photon final state 
and treating hard and soft contribntions separately. Rather than assnming that the decay 
occnrs via annihilation of two at-rest qnarks, the anthors smear the annihilation over a size 
of order 1/m, with a corresponding non-zero velocity. Althongh their calcnlation resnlts in 
some softening of the photon spectrnm relative to the lowest-order QCD prediction, it is 
nnable to entirely acconnt for the softening observed in data, leading to the conclnsion that 
final-state glnon interactions are important, particnlarly near the photon endpoint. 

Fleming and Leibovich^^ considered the photon spectrnm in three distinct momentnm 
regions. At low momentnm {x^ < 0.3), final-state radiation effects dominate. In the inter¬ 
mediate momentnm regime (0.3< x^ <0.7), they applied the operator prodnct expansion 
(OPE) to the direct photon spectrnm of T decay, with power-connting rules prescribed by 
non-relativistic QCD (NRQCD), and retained only the lowest-order color singlet terms in 
v/c. In the highest-momentum regime {x^ > 0.7), a soft-collinear effective theory (SCET) 
for the light degrees of freedom combined with non-relativistic QCD for the heavy degrees 
of freedom was used to obtain a prediction for the photon spectrum which qualitatively 
described the essential features of the CLEO-II data, despite peaking at a higher value 
of x^ than data. The same approach was later applied by Fleming to decays of the type 
—*• JZ-^-l-X, given the similarity to ^ T(1S)—>• 7 -I-X, and including the color-octet 
contributions to J/^jj production [l^. 

Very recently, Garcia and Soto (gs[I 3 ) have also produced a parameterization of the 
expected photon momentum spectrum in the T system. Following Fleming and Leibovich, 
they also remedy the inability of non-relativistic QCD (NRQCD) to model the endpoint 
region by combining NRQCD with Soft-Collinear Effective Theory, which allows calculation 
of the spectrum of the collinear gluons resulting as ^ 1. They make their own calculation 


^ It is important to note here that all these measurements assumed that the Koller-Walsh angular distribu¬ 
tion was still applicable to the phenomenological Field model. 

^ We emphasize here that the CLEO-II analysis, in presenting the background-subtracted direct photon 
spectrum, showed only statistical errors, whereas the systematic errors in the region x-y < 0.4 are con¬ 
siderably larger than those statistical errors. This is a point that was not made strongly enough in the 
past, encouraging various theoretical fragmentation models to be tested against those direct photon data 
at low x-y values. 
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of the octet contributions (in both S- and P- partial waves) to the overall rate, obtaining 
a spectral shape prediction similar to Fleming and Leibovich (claimed to be reliable in the 
interval 0.65< <0.92^3) after adding color-octet, color-singlet, and fragmentation con¬ 

tributions. For <0.65, Garcia-Soto consider the fragmentation contribution “signihcant” 
compared to the direct photon spectrum. For x^ >0.92, the calculation becomes less reli¬ 
able; in this high-momentum regime, the possibility of two-body decays: T — gg'y —>■ dfy, 
with X some resonant hadronic state, will also lead to distortions of the expected spectrum. 
Contributions from such possible two-body decays may also result in a slight underestimate 
of the extracted value of as- Garcia and Soto have also pointed out the possibility of different 
calculational regimes for T(1S)^ 7 X, compared to T(2S)^ 7 -|- X and T(3S)—>• 7 -|- X, 
given the difference in the principal quantum numbers, and therefore the average radial 
interquark separation. Since the Field model is based on simple gluon-gluon fragmentation 
phase space arguments, it does not distinguish between direct photons from any of the three 
T resonances. 


Overview of Analysis 

The analysis, in general terms, proceeds as follows. After selecting a high-quality sample 
of e’''e“ annihilations into hadrons, we plot the inclusive isolated photon spectrum in data 
taken at both on-resonance and off-resonance energies. A direct subtraction of the off- 
resonance contribution isolates the photon spectrum due to T decays. The background from 
decays of neutral hadrons into photons (tt® — 77 , rj 77 , rj' —>• p 7 , and oj —> 7 r° 7 ) produced 
in T decays to ggg, gg'y, or qq is removed statistically using a Monte Garlo generator 
developed specihcally for this purpose, and based on the assumption that the kinematics 
of charged and neutral hadron production can be related through isospin conservation. For 
each charged pion identihed in the data, we simulate a two-body decay of one of the neutral 
hadrons enumerated above. The measured four-momentum of that charged pion is then used 
to boost the daughter photons into the lab frame. After correcting for efficiency, and scaling 
by the expected rate of neutral hadron production relative to charged pion production (for 
tt'^’s, the simple isospin assumption would be X( 7 r‘’)/X( 7 r^) ~l/ 2 ; for the other neutral 
hadrons, we use ratios relative to charged pions as derived in our previous analysisjlH) and 
the appropriate branching fractions, a background “pseudo-photon” spectrum is created. 
After subtracting all backgrounds, the remaining photon spectrum is interpreted as the 
direct photon spectrum, which must then be extrapolated into low-photon momentum and 
high cos 6z regions (for which the backgrounds are prohibitively large) in order to determine 
an estimate of the full production rate. In this analysis, we employ the models by Field 
and Garcia-Soto for integration purposes, given their acceptable match in spectral shape 
to previous data. Although no predictions exist for direct photon decays of the T(2S) and 
T(3S) resonances, we nevertheless use these same models to determine total direct photon 
decay rates in the case of these higher resonances. A comparison of the shapes of these 
models is shown in Figure ^ 


Data Sets and Event Criteria 

The GLEG III detector is a general purpose solenoidal magnet spectrometer and calorime¬ 
ter. Elements of the detector, as well as performance characteristics, are described in detail 
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FIG. 1: Comparison of the direct photon spectral shapes for the two theoretical models used in 
this analysis. 


elsewhere [l^ . 11^ . 1^ . For photons in the central “barrel” region of the cesinm iodide (Csl) 
electromagnetic calorimeter, at energies greater than 2 GeV, the energy resolution is given 
by 

= + ( 2 ) 

where E is the shower energy in GeV. At 100 MeV, the calorimetric performance is about 
20% poorer than indicated by this expression due to the material in front of the calorimeter 
itself. The tracking system, Ring Imaging Gerenkov Detector (RIGH) particle identihcation 
system, and electromagnetic calorimeter are all contained within a 1 Tesla superconducting 
coil. 


Event Selection 

The data used in this analysis were collected on the T(1S) resonance, center-of-mass 
energy Eqm = 9-46 GeV, the T(2S) resonance, center-of-mass energy Eqm = 10.02 GeV, 
and the T(3S) resonance, center-of-mass energy Eqm = 10.36 GeV. In order to check our 
background estimates, we used continuum data collected just below the T(1S) resonance, 
center-of-mass energy 9.431 GeV < Eqm < 9.434 GeV, below the T(2S) resonance, center- 
of-mass energy 9.996 GeV < Ecm < 10.004 GeV, below the T(3S) resonance, center-of-mass 
energy 10.329 GeV < Eqm < 10.331 GeV and below the T(4S) resonance, center-of-mass 
energy 10.41 GeV < Eqm < 10.57. 

To obtain a clean sample of hadronic events, we selected those events that had a minimum 
of four high-quality charged tracks (to suppress contamination from QED events), a total 
visible energy greater than 15% of the total center-of-mass energy (to reduce contamination 
from two-photon events and beam-gas interactions), and an event vertex position consistent 
with the nominal collision point to within ±5 cm along the axis {z) and ±2 cm 
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in the transverse (r — 0) plane. We additionally veto events with a well-defined electron or 
muon, or consistent with a rr “1-prong vs. 3-prong” charged-track topology. Our full data 
sample is summarized in Table I. 


DataSet 

Resonance 

C (pb-^) 

HadEvts 

Wr (nb) 

EvtSel (raw) 

IS-A 

T(1S) 

6.351 

128019 

20.16 

226746 

IS-B 

T(1S) 

633.399 

12803279 

20.21 

15720815 

IS-C 

T(1S) 

424.668 

8742850 

20.59 

10553140 

2S-A 

T(2S) 

450.907 

4165745 

9.24 

6561803 

2S-B 

T(2S) 

6.133 

55834 

9.10 

76208 

2S-C 

T(2S) 

199.665 

1839390 

9.21 

2748240 

2S-D 

T(2S) 

248.473 

2299910 

9.26 

2914640 

2S-E 

T(2S) 

283.890 

2629250 

9.26 

3473320 

3S-A 

T(3S) 

382.902 

2482170 

6.52 

3887570 

3S-B 

T(3S) 

607.122 

3948690 

6.50 

5736980 

3S-C 

T(3S) 

180.758 

1168980 

6.47 

2108220 

IS-CO-A 

< T(1S) 

141.808 

485790 

3.43 

619060 

IS-CO-B 

< T(1S) 

46.600 

159959 

3.43 

260599 

2S-CO-A 

< T(2S) 

153.367 

472071 

3.08 

624505 

2S-CO-B 

< T(2S) 

106.409 

326371 

3.07 

465939 

2S-CO-C 

< T(2S) 

32.153 

99377 

3.09 

138898 

2S-CO-D 

< T(2S) 

59.783 

183897 

3.08 

256185 

2S-CO-E 

< T(2S) 

44.635 

137083 

3.07 

191205 

3S-CO-A 

< T(3S) 

46.906 

135069 

2.88 

193749 

3S-CO-B 

< T(3S) 

78.947 

226700 

2.87 

321169 

3S-CO-C 

< T(3S) 

32.064 

91997 

2.87 

130021 

4S-CO-A 

< T(4S) 

215.604 

594662 

2.76 

847875 

4S-CO-B 

< T(4S) 

558.442 

1536020 

2.75 

2189720 

4S-CO-C 

< T(4S) 

270.896 

753418 

2.78 

1073410 

4S-CO-D 

< T(4S) 

656.261 

1815920 

2.77 

2587650 

4S-CO-E 

< T(4S) 

238.903 

660883 

2.77 

941162 

4S-CO-F 

< T(4S) 

338.620 

938454 

2.77 

1337990 


TABLE I: Summary of data used in analysis. Different running periods are designated by capital 
roman letters. For each data set, we track the number of photons per unit luminosity, as well as 
the total number of observed hadronic events per unit luminosity; consistency of our results across 
datasets is later used as part of our systematic error assessment. EvtSel denotes events analyzed, 
HadEvts denotes the total number of events in each sample identified as hadronic by our event 
selection requirements, and a is the corresponding observed hadronic cross-section for each data 
sample. 


Determination of Ngg^ 

To obtain Ngg^, we had to determine the number of direct photon events. Then, with 
the number of three-gluon events, we can extract the ratio R^. For Ngg^, only photons 
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from the barrel region (| cos 6 * 2 ! < 0.7) were considered. Photon candidates were required to 
be well-separated from charged tracks and other photon candidates, with a lateral shower 
shape consistent with that expected from a true photon. Photons produced in the decay 
of a highly energetic vr® would sometimes produce overlapping showers in the calorimeter, 
creating a so-called ‘merged’ vr*^. Two selection requirements were imposed to remove this 
background. First, any two photons which both have energies greater than 50 MeV and also 
have an opening angle such that >0.975 are removed from candidacy as direct 

photons. Second, an effective invariant mass was determined from the energy distribution 
within a single electromagnetic shower. Showers with effective invariant masses consistent 
with those from merged vr^’s were also rejected. After all photon and event selection re¬ 
quirements, the momentum-dependent direct-photon hnding efficiency is shown in Figure |2l 
as calculated from a large-statistics sample of photon showers simulated with the standard, 
GEANT-based CLEO III detector simulation. We note that, since the minimum charged 
multiplicity requirement dominates the efficiency near the upper endpoint, the T(2S)—> gg^y 
and T(3S)—>■ ggj direct-photon hnding efficiencies e{x^) are higher than those shown for 
the T(1S), given their higher initial center-of-mass energies. Our hnal branching fraction 
calculations explicitly correct for this photon momentum dependence. 



FIG. 2: Efficiency for an event containing a fiducially contained direct photon to pass both event 
selection and shower selection requirements for T(1S)—> 'jgg using our default photon selection 
requirements and our default charged multiplicity requirement (>4 charged tracks observed in a 
candidate event; solid); T(1S)—> jgg showing the efficiency if the multiplicity requirement was 
relaxed to >2 charged tracks (dashed); T(1S) direct-photon daughters, for T(1S) produced in 
T(2S) dipion decays (dotted, and discussed later in this document). Efficiencies are derived from 
full GEANT-based GLEO III detector simulations. 

Using GEANT-based GLEO III detector simulations, we have compared the shower- 
reconstruction efficiency (not imposing event selection requirements) for direct photons with 
the shower-reconstruction efficiency for well-separated photons produced in the decay 7 r° —> 
77 ; these efficiencies are observed to agree to within 3% over the momentum region of 
interest (Fig. O]). The photon spectrum inferred from the observed charged pion spectrum is 
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“multiplied” by the dashed line in this hgnre to estimate the backgronnd photons expected 
from the decay of nentral pions prodnced in glnon and qnark fragmentation. 



FIG. 3: Comparison of shower-reconstruction efficiency for direct photons, compared to well- 
separated photons resulting from decays, based on full GEANT-based CLEG III simulations. 
The former efficiency is used to determine final direct photon signal branching fractions in the 
data itself; the latter is used to determine the fraction of generated pseudo-photons which are 
expected to contribute to the background showers observed in data. The difference between the 
two is attributed to the typically greater isolation of direct photons. 

The dominant backgronnds to the direct photon measnrement are of two types: initial 
state radiation {x.y > 0.65) and the overwhelming nnmber of backgronnd photons primarily 
from asymmetric 7 r° decays <0.65), that result in two, spatially well-separated daughter 
photons which elude the suppression described above. If our Monte Carlo generator 

were sufficiently accurate, of course, we could use the GEANT-based CLEO III Monte 
Carlo simulation itself to directly generate the expected background to the direct photon 
signal, including all background sources. We have compared this GEANT-derived photon 
spectrum (based on the JETSET 7.4 event generator) with data for continuum events at 
Ecu = 10.55 GeV. We observe fair, but not excellent agreement between the two, motivating 
a data-driven estimate of the background to the direct photon signal. We use GEANT to 
model the response of the calorimeter to photon showers (Fig. |2I), but use the data itself 
as an event generator of three-gluon decays, in place of JETSET. To model the production 
of 7 r° daughter photons, we took advantage of the similar kinematic distributions expected 
between charged and neutral pions, as dictated by isospin invariance. Although isospin 
conservation will break down at low center-of-mass energies (where, e.g., the neutral vs. 
charged pion mass differences and contributions from weak decays may become important), 
at the high-energy end of the spectrum (provided there is sufficient phase space), we expect 
isospin conservation to be reliable, so that there should be half as many neutral pions as 
charged pions. We stress here that this is true for three-gluon decays, XbJ decays, 1 = 0 
continuum qq events, 7 = 0 e+e“ —> T —> 7 * —> gg events, etc. 
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There are, nevertheless, both ‘physics’ as well as detector biases which comprise cor¬ 
rections to our isospin assumption, as follows. For continuum production of hadrons via 
7 * —iuu + dd), the ratio of / = 1// = 0 production is expected to be 9:1. Particles with 
7=0 (a;, /o, etc.) should decay in accordance with our naive assumption that 7 r°/ 7 r^=l/ 2 . 
For sufficiently high-multiplicity decays, such that all p states are populated evenly, we 
again expect 7 r‘^/ 7 r^=l/ 2 ; very close to the threshold turn-on, phase space effects will favor 

production, in which case <1/2. Our explicit subtraction of the photon spectrum 

obtained on the continuum will remove any such biases from e’''e“ —*• 7 * ^ gg, leaving 
three-gluon decays as the primary background source, which are presumed to obey isospin 
conservation. 

Acceptance-related biases, which will affect both continuum as well as resonance decays, 
include: a) slight inefficiencies in our charged tt^ identification and tracking, b) charged 
kaons and protons which fake charged pions, and c) for low multiplicity events, an enhanced 
likelihood that an event with charged pions will pass our minimum charged-multiplicity 
requirement compared to an event with neutral pions. The ratio therefore deviates 

slightly from 0.5, as a function of momentum. Figure 0 shows the (JETSET-I-GEANT)- 
based neutral to charged pion production ratio for continuum gg( 7 ) events taking 

into account such selection biases; we observe agreement with the 0.5 expectation to within 
3%. For this study, we rely on JETSET 7.4 to produce the proper ratio of vr^ : tt^ at 
the generator level in hadronic T fragmentation, if not the individual spectra themselves. 
In our analysis, we use this ratio, rather than the simple isospin expectation, to generate 
pseudo- 7 r°’s using data charged pions as input. The deviation between this value and the 
simple isospin expectation is later incorporated into the overall systematic error. 



FIG. 4: 7r^/(7r+-|-7r“) ratio, as a function of charged pion momentum, including tracking efficiency, 
particle identification efficiency and event selection requirements (from GEANT-based CLEG III 
Monte Garlo simulations). The loss of efficiency at large p-j^ is largely due to the bias introduced 
by the minimum charged-particle multiplicity requirement. 

These pseudo-vr^’s are subsequently decayed according to a phase space model, and the 
resulting simulated photon spectrum is then plotted. It includes our GEANT-derived photon 


10 


















efficiency, and the correlation between daughter photon momenta and the photon emission 
direction relative to the flight direction in the lab. In addition to tt^’s, we also simulate 


r; ^ 77, a; 


7r°7, 


and rj' 


Mp , 00 , 7 ) contributions, using previous measurements of these 
backgrounds in T(1S) decays 1^. An estimate of the relative contribution of these various 
backgrounds to the observed continuum spectrum is shown in Figure El 



FIG. 5: Estimate of the momentum-dependent contribution from various background sources to 
the observed below-T(4S) inclusive photon spectrum, based on the JETSET 7.4 event generator 
plus a full CLEO III GEANT-based detector simulation. Initial state radiation contributions are 
not included. 

We have also studied the relative contribution to the inclusive spectrum from neutrons, 
anti-neutrons, and iF°’s. According to Monte Carlo simulations, the expected numbers of 
such particles per hadronic event with | cos( 6 '^)| < 0.7, and scaled momentum x > 0.25 {i.e., 
particles which could populate our signal region) are quite small. Figure El gives the yield 
per event, as a function of momentum, for and anti-neutrons to contaminate our signal 
region. Such contributions are therefore neglected in the remainder of the analysis. 

The performance of our photon-background estimator can be calibrated from data itself. 
Three cross-checks are presented below: a) comparison between the absolutely normal¬ 
ized angular distribution of our simulated pseudo-photons (“PP”) using continuum charged 
tracks as input to our pseudo-photon generator^ versus the photon spectrum measured on 
the continuum (including a Monte Carlo-estimated initial state radiation (“ISR”) contribu¬ 
tion (Figure 13) ), b) comparison of the absolute magnitude of the pseudo-photon momentum 
spectrum with continuum data (Figure 13,^ and c) comparison of the reconstructed 7 r° and 
rj mass peaks (Figure Ej) between our simulated photons and real data photons. All these 


^ Note that there are two simulations referred to in this document - “simulated” PP photons refer to the 
pseudo-photons generated using identified charged pion tracks as inputs; “Monte Carlo” refers to the full 
JETSET-l-GEANT CLEO III event-l-detector simulation. 

^ Note that the Monte Carlo ISR (MC ISR) spectrum shows an enhancement in the interval 0.7 < Xj < 0.8, 
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FIG. 6: Estimate of the momentum dependent percentage of showers produced by h^’s and 
that passed our shower selection, based on a sample of 1 million MC continuum events. There are 
no entries for x > 0.5. 

checks show acceptable agreement between simnlation and data. The nnmerical accnracy 
of onr backgronnd estimate can be assessed by comparing, for the second of these checks, 
the fractional excess remaining after the estimated psendo-photon backgronnd (+ISR) is 
snbtracted from the raw continunm data spectrnm, in the momentum interval of interest: 
(A^data - A^pseudo-photons “ A^isr)/A^ data- Integrated from a;^=0.4 to x.^=0.95, we hnd the frac¬ 
tional excesses to be -1.86%, -0.68%, 2.55% and 1.76%, using data iDelow the T(1S), T(2S), 
T(3S) and T(4S) resonances, respectively; we consider these excesses to be acceptably con¬ 
sistent with zero. 

Other systematic checks of our data (photon yield per data set, comparison between the 
below-T(lS), below-T(2S), below-T(3S) and below-T(4S) continuum photon momentum 
spectra) indicate good internal consistency of all data sets considered. 


Signal Extraction 

Two different methods were used to subtract background photons and obtain the 
T(1S)—>• gg'y^ T(2S)—> gg'y and T(3S)—> gg'y spectra. In the hrst, after explicitly sub¬ 
tracting the continuum photon spectrum from data taken on-resonance, we use the pseudo¬ 
photon spectrum to model the background due to 7r°, rj, r]' and a; decay which must be 
separated from direct photons from T decay. In the second method, we used an exponential 
parametrization of the background to estimate the non-direct photon contribution. 


compared to the lack of events in the region 0.8 < x-y. This is attributable to a) the e“'"e“ —> ccy 
threshold being crossed for x-y > 0.8, b) since the e+e“ —> qq cross-section aqq ^ 1/s ~ 1/a:^, there is an 
enhancement in hadronic final-state production as the energy of the radiated ISR photon approaches the 
beam energy. 
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FIG. 7: Angular distribution of the inclusive photons from continuum data compared with our 
pseudo-photon estimate, based on isospin invariance, for showers with > 0.45, and including a 
Monte Carlo-based estimate of the ISR background. The normalization is absolute. 


3030805-007 



FIG. 8: Comparison of spectra, obtained using below-T(4S) data, with the sum of ISR Monte 
Carlo simulations plus a pseudo-photon spectrum obtained using identified data charged pions as 
input. 


Figure ^1 shows the inclusive T(1S) photon distribution with the different estimated 
background contributions (continuum photons from all sources and T decays of neutral 
hadrons into photons) overlaid. After subtracting these sources, what remained of the 
inclusive T(IS) spectrum was identihed as the direct photon spectrum, T(IS)— gg'-^. Figures 
HH and ini show the corresponding plots for the T(2S) and T(3S) data, and also indicate 
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FIG. 9: The and rj yields for data (dashed) and simulated photons (solid). The yields agree at 
the 2%-3% level. 


the magnitude of the cascade subtraction due to transitions of the type T(2S)—> T(1S)+X, 
T(1S)—> '^gg. We use currently tabulated values for T(2S)—>• T(1S)+X to determine the 
magnitude of this correction. Monte Carlo simulations of the primary cascade processes, 
including T(2S)—^ 7r7rT(lS) (using a Yan 
and T(2S)^ Xbl, Xb 
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distribution for the dipion mass distribution) 
7T(1S) are used to adjust the shape of our measured T(1S)^ ■ygg 
direct photon spectrum to that expected for the cascade subtraction in order to account for 
the shifted kinematic endpoint and Doppler smearing of the daughter T(1S) direct photon 
spectrum. We assume that the daughter T(1S) retains the polarization of the parent T(2S); 
the direct photon angular distribution is then the same as for direct production and decay 
of the T(1S) resonance. 


Parametric estimate of background 

Observing that the photon spectrum seems to describe an exponential outside the signal 
region, we attempted to check our pseudo-photon and continuum-subtracted yields against 
the signal photon yields obtained when we simply fit the background to an exponential in 
the momentum region below the signal region (comparing the results obtained from fitting 
0.2 < < 0.3 to those obtained using 0.3 < < 0.4) and then extrapolated to the region 

0.4 < x^. Figure UHl shows that this procedure satisfactorily reproduces continuum data 
below the T(4S) resonance, verifying that it may be used to generate a rough estimate of 
the backgrounds. 
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FIG. 10; Photon energy spectrum (minimal cuts) for data taken at the T(1S) resonance energy, 
with continuum contribution (using CLEO data taken off the resonance, including ISR) and T 
non-direct simulated pseudo-photons (“pp”) resulting from decays of neutral hadrons overlaid. 


3030805-004 



FIG. 11: T(2S) photon energy spectrum (minimal cuts), with non-direct pseudo-photons, contin¬ 
uum background photons, and the cascade contribution from T(1S) decays overlaid. 

Model Fits 

We estimate by extrapolating the background-subtracted photon spectrum down to 
using a model to prescribe the spectral shape at low photon momentum. Since the 
CLEO calorimeter has hnite resolution, and since the photon-hnding efficiency is momentum- 
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FIG. 12: T(3S) photon energy spectrum (minimal cuts), with non-direct pseudo-photons, contin¬ 
uum background photons, and the cascade contributions from T(2S) and T(1S) decays overlaid. 



Xy 


FIG. 13: Subtraction of backgrounds using an exponential (below-4S continuum data), with float¬ 
ing normalization to estimate the non-direct photon spectrum. The exponential, plus initial state 
radiation, gives a fair match to the observed spectrum, although the background is clearly under¬ 
estimated in the intermediate region of the momentum spectrum. 


dependent, two procedures may be used to compare with models. Either a migration-matrix 
can be determined from Monte Carlo simulations to estimate the bin-to-bin smearing, with 
a matrix-unfolding technique used to compare with prediction, or the model can first be 
efficiency-attenuated (as a function of momentum) and then smeared by the experimental 
resolution to compare with data. We have followed the latter procedure, floating only the 
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normalization of the efficiency-attenuated, resolution-smeared model, in this analysis. To 
determine the percentage of direct photons within our hducial acceptance, we used the QCD 
predictions of Koller and Walsh for the direct-photon energy and angular distributions!^. 
Our large statistics sample allows (for the hrst time) a check of the Koller-Walsh predic¬ 
tion. Figure ITTI shows that the angular distribution of our data, after taking into account 


3031005-026 



cos 0 


FIG. 14: Photon angular distribution for background-subtracted direct photon data (histogram) 
vs. Koller-Walsh prediction, modified for the experimental efficiency as a function of and cos 6z- 

acceptance effects, agrees adequately with the Koller-Walsh prediction. 

Figures ITHl inn and El show the hts of the direct photon energy spectrum to the Garcia- 
Soto direct photon model. The hts are performed over the interval claimed to be relatively 
free of either endpoint effects or fragmentation backgrounds (0.65< <0.92), then extrap¬ 

olated under these backgrounds into the unht region using only the direct photon component 
of their spectral model. Field prescribes no such cut-offs, so we have ht that model over the 
larger kinematic range 0.4< <0.95. To probe htting systematics, we have performed two 

hts. In the hrst, we perform a simple minimization of the background-subtracted data 
to the Garcia-Soto spectrum. In the second, we have normalized the area of the theoretical 
spectrum to the area of the background-subtracted data in the interval of interest. The two 
methods yield nearly identical results. 

We note, in some cases, an excess of photons in data as a;.y —> 1. Further examination of 
these events indicate that they are dominated by e’''e“ —»■ 77r+7r“7r+7r“. 

Figures HSl HE and 1201 show hts obtained using a simple exponential parametrization of 
the background, with no pseudo-photon generation. 


T(2S)^ 7r+7r-T(lS); T(1S)^ 795 

Our large sample of T(2S) decays and the substantial T(2S)—>• 7r+7r“T(lS) branching 
fraction (~0.19) ahord an opportunity to measure a ‘tagged’ T(1S) direct photon spectrum 
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FIG. 15; Fit to background-subtracted T(1S) data, using explicit continuum data subtraction. 
Garcia-Soto model is used for spectral shape (modified for efficiency and experimental resolution), 
either using a fit the region where the direct photon contribution dominates, or normalizing 
the model to the experimental data in the same interval, as shown. The two fits very nearly overlay 
with each other. 

which circumvents all continuum backgrounds. In a given event taken at the T(2S) center- 
of-mass energy, we calculate the mass recoiling against all oppositely-signed charged pion 
pairs fFigure I71|l . In each bin of recoil mass, we plot the spectrum of all high-energy photons 
in that event. A sideband subtraction around the T(2S)—> 7r’''7r“T(lS) recoil mass signal, 
at rUfecoil rsj T(1S) results in a tagged T(1S) direct photon spectrum. Spectral shape and 
R.y values obtained this way are consistent with our other estimates. 


Determination of N, 


999 


To determine the number of three-gluon events Nggg from the number of observed T(1S) 
hadronic events we first subtracted the number of continuum events at the T(1S) en¬ 
ergy based on the observed number of below-T(lS) continuum events 


f=9.43 GeV 


^T(IS) ^ ^_Bcm=9.43 GeV . '^T(lS)-Scont _ 
cont cont ' n 

•^cont'^lS 


coni 


( 3 ) 


where the factor arises from the number of events N = C ■ and s = Eqj^. From 

the observed number of hadronic events collected at the resonance and knowing the 

branching fractions and efficiencies for T(1S) ^ qq, T(1S) ^ 99 ’^-, and T(1S) —> 999 ^ the 
number of T(1S)—>• 999 events can be inferred. For TflS) —>• gg, e.g., we use the averaged 
T(1S)—> branching fraction = 0.0248 [13, and Rr{is) = c’'(e+e“ —T(1S)—>• 

qq)/a{e^e~ T(1S)— = 3.51 jl^, 

.T(1S) 


Nqq — R 


T(1S) ■ Bg^eqq 


N, 


had 


(1 


( 4 ) 
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FIG. 16: Fit to background-subtracted T(2S) data, using explicit continuum data subtraction and 
explicit subtraction of T(1S) cascade contributions. Direct spectrum fit using Garcia-Soto model. 



FIG. 17: Fit to background-subtracted T(3S) data, using explicit continuum data subtraction and 
explicit subtraction of T(1S) and T(2S) cascade contributions. 


with Shad the hadronic event reconstruction efficiency, averaged over all hadronic modes, and 
eqq specihcally the efficiency for reconstructing an T(1S) —(gg) event. 

Using B{T{2S)^ T(1S) + X) = (32±l)%, i3(T(3S)^ T(2S) + X) = (10.6±0.8)%, and 
i3(T(3S)—> T(3S) -I- X) = (12.1±0.5)%, three-gluon decay fractions for the three resonances 
fggg=81.3 0.5%, 39 it 1%, aud 38 ± 1%, respectively, we obtained a value for Nggg for each 

of the resonances, based on our measured values for Ngg^. More details on this subtraction 
are presented in Table II. 
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Event Type 

Efficiency (e) 

Tr(lS)^ ggg 

0.953 ± 0.003 

Tr(lS)^ ggX 

0.751 ± 0.007 

T{lS)^qq 

0.871 ± 0.005 

T(2S)^ ggg 

0.956 ± 0.003 

T(2S)^ ggX 

0.776 ± 0.007 

T(2S)^ qq 

0.882 ± 0.005 

T{2S)^T{IS)+X ^ggg 

0.956 ± 0.003 

T(2S)^ T(IS)+X ^ ggx 

0.778 ± 0.007 

T(2S)^ T(1S)+A ^ qq 

0.891 ± 0.005 

T(2S)^ Xbj{lV)^gg (J = 0,l,2) 

0.933 ± 0.004 

T(3S)^ ggg 

0.955 ± 0.003 

T(3S)^ ggX 

0.765 ± 0.007 

T(3S)^ qq 

0.881 ± 0.005 

T(3S)^T(2S)+X^5ra 

0.958 ± 0.003 

T(3S)^ T(2S)+X ^ ggx 

0.765 ± 0.007 

T(3S)^ T(2S)+A ^ qq 

0.877 ± 0.005 

Ti3S)^TilS)+X ^ ggg 

0.961 ± 0.003 

T(3S)^ T(1S)+X ^ ggj 

0.789 ± 0.006 

T(3S)^ T(1S)+A ^ qq 

0.90 ± 0.07 

T(3S)^ Xbjm^gg {j = 0,1,2) 

0.819 ± 0.006 

T(3S)^ Xbj{2P)^gg (J = 0,I,2) 

0.929 ± 0.004 

T Resonance 

iVtotai(T(nS)) (10^) 

T(1S) 

21.0 ± 0.06 

T(2S) 

8.4 ± 0.04 

T(3S) 

5.2 ± 0.06 

Traction 

/ 

/(T(1S)^ ggg) 

0.813 ± 0.005 

/(T(2S)^ ggg) 

0.39 ± 0.01 

/(T(3S)^ ggg) 

0.38 ± 0.01 

/(T(2S)^ T(1S)+X) 

0.32 ± 0.01 

/(T(3S)^ T(2S)+X) 

0.106 ± 0.008 

/(T(3S)^ T(1S)+X) 

0.121 ± 0.005 


TABLE II: Efficiencies for the reconstruction of the various types of events considered in this 
analysis, the total number of calculated T(1S), T(2S) and T(3S) events, and the fractions of these 
totals which were used to obtain the T —> ggg denominator in our measurements of (and to 
scale the direct photon cascade spectra in our T(2S) and T(3S) subtractions). These fractions 
were obtained from the Particle Data Group [i^ and include recent CLEO Xb.j measurements [ 2 ^. 
The presented errors on the efficiencies are statistical only. Note that, although there are more 
possible decay paths from the T(3S) than the T(2S), the ggg fractions are comparable owing to 
the significantly larger T(2S)^ T(lS)7r7r branching fraction. 
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FIG. 18: Subtraction of backgrounds using an exponential (T(1S) data), with floating normaliza¬ 
tion to estimate the non-direct photon spectrum. Direct spectrum fit using Field model. 
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FIG. 19: Subtraction of backgrounds using an exponential (T(2S) data), with floating normaliza¬ 
tion to estimate the non-direct photon spectrum. Direct spectrum fit using Field. 


Results 

With values of Ngg^ and Nggg, the ratio can be determined. Table III presents our 
numerical results for the extracted branching fractions. We note that, in general, the reduced 
values for the hts tend to be rather high. Structure in the spectrum due to, e.g., two- 
body radiative decays, may result in such a poor £t and is currently being investigated. The 
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FIG. 20: Subtraction of backgrounds using an exponential (T(3S) data), with floating normaliza¬ 
tion to estimate the non-direct photon spectrum. Direct spectrum fit using Field. 



9.42 9.44 9.46 9.48 9.50 

^recoil 

FIG. 21; Mass recoiling against oppositely signed charged pion pairs, T(2S) data. 

normalization-by-area fits probe the extent to which the model hts may be disproportionately 
weighted by a small number of points. 

Systematic errors 

We identify and estimate systematic errors as follows: 
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A ^ ggr, X= 

Background 

Pield R^ / x^/d.o.f. 

GS R^ / xVd.o.f. 

T(1S) 

T(1S) 

T(1S) 

T(2S) ^ 7r7rT(lS)-tagged 
T(2S) 

T(2S) 

T(2S) 

T(3S) 

T(3S) 

T(3S) 

Exponential 
PP (MC ISR) 
PP (CO ISR) 
PP (no ISR) 
Exponential 
PP (MC ISR) 
PP (CO ISR) 
Exponential 
PP (MC ISR) 
PP (CO ISR) 

(2.94 ±0.02)% / 115.1/67-1 
(2.81 ± 0.01)% / 293/74-1 
(2.93 ± 0.01)% / 125/67-1 
(2.9 ±0.3)% / 118/58-1 
(3.7 ± 0.7)% / 542/105-1 
(3.42 ± 0.05)% / 316/67-1 
(3.58 ±0.05)% 145/67-1 
(3.4 ±0.4)% 210/105-1 
(2.91 ± 0.07)% / 263/67-1 
(2.8 ±0.1)% / 72/67-1 

(2.39 ± 0.03)% / 132.4/67-1 
(2.48 ± 0.01)% / 694/74-1 
(2.45 ±0.01)% / 116/37-1 
(2.5 ± 0.3)% / 132/58-1 
(3.4 ± 0.4)% / 773/105-1 
(3.01 ± 0.04)% 426/67-1 
(2.77 ± 0.05)% 87/37-1 
(3.1 ±0.1)% 251/105-1 
(2.55 ±0.06)% 331/67-1 
(2.1 ±0.1)% / 36/37-1 


TABLE III: Summary of Measurements. “PP” denotes Pseudo-Photon background, “MC ISR” 
implies that Monte Carlo simulations of initial state radiation were used to subtract the ISR 
background. These numbers are provided for comparison only and are not used in final averaging, 
etc. “CO ISR” implies that data ISR was subtracted directly using below-resonance data. Values 
obtained using an exponential parametrization have had systematic errors (reflecting sensitivity to 
region chosen for scale normalization of the exponential outside the peak region; for this estimate, 
branching fractions were compared using the regions 0.2 < x-f < 0.3 or 0.3 < < 0.4 to set 

the scale of the exponential and extrapolate under the signal in the higher-a;-,, region) added in 
quadrature with the statistical error. All other errors are statistical only. Note that T(2S) 
values have been corrected for T(1S)^ gg'j contamination; T(3S) values have been corrected 
for both T(2S)^ gg^, and also T(1S)^ ggj contamination. 

1. For the T(1S), the uncertainty in Nggg is based on the CLEO estimated three-gluon 
event-finding efficiency uncertainty. For the T(2S) and T(3S) decays, the uncertainty 
in Nggg ulso folds in uncertainties in the tabulated radiative and hadronic transition 
decay rates from the parent T’s, which are necessary for determining Nggg as well as 
the magnitude of the cascade subtractions. The cascade subtraction errors include 
statistical (la) uncertainties in the various decay modes of the T resonances. 

2. Background normalization and background shape uncertainty are evaluated redun¬ 
dantly as follows: 

(a) We determine the branching fractions with and without an explicit 7r° veto on 
the background. 

(b) We measure the internal consistency of our results using different sub-samples of 
our T(1S), T(2S) and T(3S) samples. 

(c) Bias in background subtraction can also be estimated using Monte Carlo simula¬ 
tions. We treat the simulation as we do data, and generate pseudo-photons based 
on the Monte Carlo identihed charged pion tracks. After subtracting the pseudo¬ 
photon spectrum from the full Monte Carlo photon spectrum, we can compare our 
pseudo-photon and ISR-subtracted spectrum with the known spectrum that was 
generated as input to the Monte Carlo detector simulation. For the T(1S), T(2S), 
and T(3S), we observe fractional deviations of -|-5.7%, -3.4%, and -2.1% between 
the input spectrum and the pseudo-photon background-subtracted spectrum. To 
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the extent that the initial state radiation estimate and the photon-hnding effi¬ 
ciency are obtained from the same Monte Carlo simulations, this procedure is 
largely a check of our generation of the pseudo-photon background and the cor¬ 
relation of decay angle with efficiency. 

(d) We extract the direct photon branching fractions using a flat isospin ratio 

of 0.5, compared to the tt® : ratio based on Monte Carlo simulations, including 

all our event selection and charged tracking and charged particle identihcation 
systematics (~0.53, Figure E]). 

(e) Our uncertainty in the on-resonance vs. off-resonance luminosity scaling, which 
determines the magnitude of the continuum subtraction is assessed as < 1 %, ab¬ 
solute. 

3. Model dependence of the extracted total decay rate is estimated by: i) determining 
the variation between hts (for a given model) performed using a minimization 
prescription, or a simple normalization of area of theoretical spectrum to data, and 
also by ii) comparing the results obtained from hts to the Field model with results 
obtained from the hts to the Garcia-Soto model. (We currently assume the Koller- 
Walsh prescription for the angular distribution is correct, and assign no systematic 
error for a possible corresponding uncertainty.) The irreducible model-dependence 
error (the difference between branching fractions obtained with the Field model vs. 
Garcia-Soto) is presented as the last error in our quoted branching fraction. 

There is currently no theoretical consensus on either the shape or magnitude of the 
fragmentation photon background to the direct photon spectrum. We assign no explicit 
systematic error to the uncertainty in this component and presume this to be already probed 
by the variation observed between models, the consistency we observe with results obtained 
from an exponential ht to the background, and the consistency observed in htting over 
different intervals of the background-subtracted spectrum. Given the currently tabulated 
upper limit on T(1S)—> y-l-pseudoscalar, pseudoscalar^ h'^h~ {B < 3 x 10“^) and the 
small branching fractions measured for other two-body exclusive radiative decays (like the 
p(1440), whose dominant modes do not have two charged tracks in the hnal state), we neglect 
distortions to the direct photon yield from exclusive two-body decays T —>■ 7 -|- df. 

Table IV summarizes the systematic errors studied in this analysis and their estimated 
effect on R^. 


Comparison with previous analyses 

Table V compares the results of this analysis with those obtained by previous experiments, 
in which the number of T(1S)—>■ gg'y events were determined using Field’s theoretical model 
only. 


Summary 

We have re-measured the T(1S)— 'ygg/T^lS)^ ggg branching fraction ratio {R^y)-, 
obtaining agreement with previous results. We also have made hrst measurements of 


24 


Source 

6Ry{lS/2S/3S) 

Difference (MC Glevel, MC analyzed) 

0.08/0.05/0.03 

Background Shape/Norm, including: 
With/Without a veto 

Isospin Assumption 

^999 

Cascade subtraction 

0.07/0.07/0.07 

0.01/0.01/0.01 

0.08/0.19/0.28 

0/0.07/0.14 

Luminosity and ^/s scaling 

0.01/0.01/0.01 

Fit systematics (norm vs. ht) 

0.01/0.05/0.01 

Total Systematic Error 

0.13/0.22/0.32 

Model Dependence (GS vs. Field) 

0.24/0.41/0.37 


TABLE IV: Systematic Errors. 


Experiment 

Ry{%) 

CLEO 1.5 (T(lS))[a 
ARGUS (T(lS))[ll]j 
Crystal Ball (T(lS))[9j 
CLEO II (T(1S))[]JJ 

2.54 ±0.18 ±0.14 

3.00 ±0.13 ±0.18 

2.7 ±0.2 ±0.4 

2.77 ±0.04 ±0.15 

CLEO III (T(1S)) 
CLEO III (T(2S)) 
CLEO III (T(3S)) 

2.70 ±0.01 ±0.13 ±0.24 

3.18 ±0.04 ±0.22 ±0.41 

2.72 ±0.06 ±0.32 ±0.37 


TABLE V: Comparison with other experiments. Errors are statistical, systematic, and model- 
dependent (Eield vs. Garcia/Soto), respectively. Central values are obtained by a direct weighted 
average (taking into account both statistical and systematic errors, and assuming errors to be 
uncorrelated, therefore yielding the most conservative estimate of the experimental precision) of 
the measurements presented in Table III. 


R.y{2S) and R.y{3S). Our results are, within errors, consistent with the naive expecta¬ 
tion that i?.y(lS)~ R^{2S)^ R-y{3S), although this equality does not hold for the recent 
CLEO measurements of for the three T resonancesAssuming an energy scale 
equal to the parent T mass, our values of R-y for T(1S)— gg'j (2.70±0.01±0.13±0.24)%, 
T(2S)^ gg'j (3.18±0.04±0.22±0.41)%, and T(3S)^ gg-f (2.72±0.06±0.32±0.37)% imply 
values of the strong coupling constant as{Mz) = (0.1114 ± 0.0002 ± 0.0029 ± 0.0053), 
(0.1026±0.0007±0.0041±0.0077) and (0.113±0.001±0.007±0.008), respectively, which are 
within errors, albeit consistently lower, compared to the current world average (Appendix 

I). 
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Appendix I) Calculation of strong coupling constant 

The decay width T — gg'y has been calculated by Lepage and Mackenzie in terms 
of the energy involved in the decay process (he., Us^Ecm), or a!s(Mx)): 


r(T^^^7) 8(7r2-9) 2 


cXg [Ml 


1 + (3.7 ±0.4) 




TT 


(A.5) 


r(T —/i+/i ) QTTaQED 

Sanghera 2^| rewrites this expression in terms of an arbitrary energy (renormalization) scale 




r(T ^ gg-^) ^ 
r(T ^ /x+p-) " 
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TT 
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TT 
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± (3.7 ±0.4) 


(A.6) 


where = (33 — 2?xj)/127r, and Uf is the number of light quark flavors which 

participate in the process (rx/ = 4 for T(1S) decays). 

Similarly, the decay width T —ggg has been calculated by Bardeen et al. and 
expressed by Lepage et al. jH, as: 


r(T ^ ggg) 

r(T ^ /X+/X-) 


10(7r^ — 9) al{Mr] 


Slnel 
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QED L 


1± 




TT 


-[(2.770 ± 0.007)/3o - (14.0 ± 0.5)] ± 


(A^ 

with do = 11 — Ch = —the charge of the b quark. Here again Sanghera |2^ 

uses the same algebraic technique to rewrite this in terms of the renormalization scale: 


r(T ^ ggg) 

r(T ^ /X+/X-) 



± Ag 


n* (h") 

TT 
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37r6o In 




BfUf ± Bi 


(A.8) 


with Ag = d , Bf = 2.770 ± 0.007, and Bi = 16.47 ± 0.58. 

»iejj “qbd ___ 

Note that the scale dependent QCD equations ()A.6fl and (IA.8fl are hnite order in a*. If 
these equations were solved to all orders, then they could in principle be used to determine 
R.y independent of the renormalization scale. But since we are dealing with calculations that 
are hnite order, the question of an appropriate scale value must be addressed. 

The renormalization scale may be dehned in terms of the center of mass energy of the 
process, /x^ = where A is some positive fraction. Since QCD does not tell us a 

priori what A should be, we must dehne the appropriate scale. One possibility would be to 
dehne /x = Ecm] that is A~i-- A number of prescriptions [ll, |28|, |^, |33 have been proposed 
in an attempt to “optimize” the scale. However, each of these prescriptions yields scale 
values which, in general, vary greatly with the experimental quantity being measured [ 2 ^ . 
We have chosen fn = lto facilitate a calculation of ag at each of the T resonance energies. 

For the T(1S) analysis, using fi = Mr(is) we hnd 


as{Mr(is)) = 0.1735 ± 0.0005 ± 0.0072 ± 0.0133; (A.9) 

for the T(2S) analysis, using fi = My{ 2 S) we hnd 

as{Mr(2S)) = 0.151 ± 0.002 ± 0.009 ± 0.017; (A.IO) 

for the T(3S) analysis, using /x = My{ 3 s) we hnd 

as{Mr(3s)) = 0.172 ± 0.003 ± 0.018 ± 0.021. (A.ll) 
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The errors are statistical, systematic and model-dependent, respectively. These calcu¬ 
lations were obtained by finding the zeroes of the ratio of Eqs. IA.6I and IA.8I given our 
measurement of ihy for each T resonance. The errors were obtained by shifting our mea¬ 
surement of by ±cr, for each of our three errors, and extracting for each relevant 
error-shifted central value. 

These results can then be extrapolated to /i = Mz using eauation lA.121 2^ | with po = My 
for each resonance. For this calculation, only the first three terms of the /3-function were 
considered lli. 


1 ^ 

log — = 

/^o 


da 


(A.12) 


lasi/io) 

This calculation for the T(1S), T(2S) and T(3S) results in the following measurements 
of «,(Mz,T(nS)): 


af^’T(is) = 0.1114 ± 0.0002 ± 0.0029 ± 0.0053, (A.13) 

af^’T(2S) = 0.1026 ± 0.0007 ± 0.0041 ± 0.0077, (A.14) 


af^’T(3S) ^ ^ Q_QQ^ ^ Q_QQ7 ^ Q_QQg_ 

Our results are systematically low compared with the average value of a 


(A. 16) 
WM 2 ) =0.119± 

0.006 obtained from many variables studied at all the LEP experiments 22|, but in better 
agreement with as{Mz) = 0.112 ± 0.003 obtained from an analysis of structure functions in 
deep inelastic scattering and with the previous CLEO measurement of as(M^,T(lS)) 
Til l. For the T(2S) and T(3S) measurements, we stress caution in interpreting these results. 


as it is (again) unclear what procedure should be used to define the renormalization scale. 

As an alternative to the extraction method outlined above, the strong coupling con¬ 
stant as can be written as a function of the QCD scale parameter Rjig, defined in the 
modified minimal subtraction scheme (MMSS)[2^. Figure |221 presents the contour plot of 

/m)- 

Similarly, the ratio of eqns. (2) and (4) above can be used to eliminate ag and provide a 
relationship between R-y, and (Figure I^Sl)- 











3031005-021 



FIG. 22; Contour plot illustrating dependence of Og on the QCD scale parameter and the 
momentum scale /^. 



FIG. 23: Contour plot illustrating relationship between R^, the QCD scale parameter and 
the momentum scale ffj,. 
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